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Abstract. The noisy Burgers equation describing for example the growth of an inter- 
face subject to noise is one of the simplest model governing an intrinsically nonequi- 
librium problem. In one dimension this equation is analyzed by means of the Martin- 
Siggia-Rose technique. In a canonical formulation the morphology and scaling behavior 
are accessed by a principle of least action in the weak noise limit. The growth morphol- 
ogy is characterized by a dilute gas of nonlinear soliton modes with gapless dispersion 
law E oc p"^/^ and a superposed gas of diffusive modes with a gap. The scaling exponents 
and a heuristic expression for the scaling function follow from a spectral representation. 



1 Introduction 

Macroscopic phenomena far from equilibrium are ubiquitous and include phe- 
nomena such as turbulence in fluids, interface and growth problems, chemical 
reactions, biological processes, and even aspects of economical and sociological 
structures. 

In recent years much of the focus of modern statistical physics and soft con- 
densed matter has shifted towards such systems. Drawing on the case of static 
and dynamic critical phenomena in and close to equilibrium, where scaling, crit- 
ical exponents, and the concept of universality have so successfully served to 
organize our understanding and to provide a variety of calculational tools, a 
similar approach has been advanced towards the much larger class of nonequi- 
librium phenomena with the purpose of elucidating scaling properties and more 
generally the morphology or pattern formation in a driven state. 

In this context the noisy Burgers equation in one dimension provides maybe 
the simplest continuum description of an open driven nonlinear system exhibiting 
both scaling and pattern formation. This equation has the form (Forster et al. 
(1976), Forster et al. (1977)) 

^ = vV'^u + XuVu + V77 , (1) 

and was in the noiseless case for 77 = originally proposed by Burgers (Burg- 
ers (1974)) in order to model turbulence in fluids; we note the similarity with 
the Navier-Stokes equation for A = — 1. Equation (1) has the form of a con- 
served nonlinear Langevin equation, du/dt = — Vj, with fluctuating current 
j = — j/Vu — (A/2)u^ -|- rj. The linear diffusive damping term vW'^u is charac- 
terized by the surface tension or viscosity v. The nonlinear convective or mode 
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coupling term AmVu is controlled by A. Finally, the equation is driven by the 
fluctuating "white noise" i] for which we assume a Gaussian amplitude distribu- 
tion 



P{r]) = exp 



~ j dxdtr]{xtf 



(2) 



and short-range correlations in space according to the correlation function 

{r]{xt)r]{x't')) = AS{x - x')S{t - t') , (3) 



characterized by the noise strength parameter A. 

In the context of modeling a growing interface the Burgers equation gov- 
erns the local slope u = Vh of a height field h (in the Monge representation) 
characterized by the much studied Kardar-Parisi-Zhang equation (Kardar et 
al.(1986), Medina et al. (1989)) 

^ = uV'h+^{Vhf+V . (4) 

In this case, which we shall adhere to in the following, is a diffusion coefficient 
or viscosity, A a nonlinear lateral growth parameter, and r] represents noise in the 
drive or the environments. In Fig. 1 wo have sketched the growth morphology in 
terms of the height and slope fields for a typical growing interface. 

The substantial conceptual problems encountered in nonequilibrium 
physics are in many ways embodied in the Burgers-KPZ equations (1) and 
(4) describing the self-afSne growth of an interface subject to annealed noise 
arising from fluctuations in the drive or in the environment. Interestingly, the 
Burgers-KPZ equations are also encountered in a variety of other problems such 
as randomly stirred fluids, dissipative transport in a driven lattice gas, the prop- 
agation of flame fronts, the sine-Gordon equation, and magnetic flux lines in 
superconductors. Furthermore, by means of the Cole- Hopf transformation the 
Burgers-KPZ equations are also related to the problem of a directed polymer or 
a quantum particle in a random medium and thus to the theory of spin glasses; 
see e.g. (Halpin-Healy and Zhang (1995)). 

One issue which has been addressed is the scaling properties of the noisy 
Burgers equation. In Fig. 2 we have depicted how the width of the fluctuating in- 
terface (see Fig. 1) after a transient lapse of time, characterized by the a crossover 
time tcrossover Scaling with the system size L according to tcrossover oc L^, where 
z is the dynamic exponent, saturates to a value wq oc also depending on L 
and characterized by the roughness exponent C. 

This dynamical scaling hypothesis, which is substantiated by numerical mod- 
eling and renormalization group considerations, is embodied in the following 
relationship for the width: 



wit,L) = L'^Fi{t/L') 



(5) 
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Fig. 1. We depict the growth morphology in terms of the height (h) and slope (m) fields 
for a growing interface. The saturated width in the stationary state is denoted by wo- 



W 




Fig. 2. We depict the interface width 'w{t) as a function of time t. In the transient 

regime t <g; tcrossover ^ L'' w grows according to t''^^. In the stationary regime attained 
for t ^ tcrossover u) Saturates to the value wo ^ . 
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In terms of the height field itself or the associated slope field the appropriate 
dynamical scaling forms are given by 

{[h{x + xo,t + to) - h{xQ, to)f} = x^^F2{tlx^) , (6) 

and 

{u{x + XQ,t + to)u{xo,to)) = x'^'^-^F3{t/x') . (7) 

Here xq and to are reference points in an infinite system in the stationary regime 

and Fi the scaling functions. The scaling issue thus amounts to a determination 
of the scaling exponents ( and z together with the scaling functions F^. 

The hydrodynamical origin of the Burgers equation, as reflected by the pres- 
ence of the mode coupling or convcctivc term XuVu, implies that the Burgers 
equation is invariant under a Galilean transformation 

X ^ X — Xucit (8) 
u ^ u + uq , (9) 

involving a shift of the slope field. Since the nonlinear coupling strength A enters 
as a structural constant in the symmetry group it is invariant under scaling. This 
property implies a dynamical scaling law 

C + ^ = 2 , (10) 

relating the roughness and dynamic exponents. 

Another interesting property of the Burgers equation specific to one dimen- 
sion is the existence of an effective static fiuctuation-dissipation theorem, in the 
sense that the stationary Fokker-Planck equation for the Btirgers equation is 
solved by a Gaussian distribution 



Pst(M) oc exp 



-^Jdxuix) 



(11) 



for the slope field independent of the nonlinear coupling strength A. This implies 
that the slope field u performs independent Gaussian fluctuations and that, 
consequently, the height field h — J udx, performs random walk in x, yielding 
the roughness exponent and from the scaling law (10) also the dynamic exponent. 
Summarizing, the universality class for the Burgers equation is characterized by 
the exponents: 

C = ^ and z=^ . (12) 

By a combination of the Galilean invariance and the effective fluctuation-dissipation 
theorem specific to one dimension the exponents in one dimension are thus ex- 
actly known. The scaling function has also been accessed numerically and by 
a mode coupling approach. In higher dimensions the scaling properties of the 
Burgers-KPZ are controversial and remain to be clarified. In this context the 
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dynamic renormalization group based on an expansion in }?Alv^ and an expan- 
sion about the critical dimension d = 2 yields limited results and in particular 
fails to access the strong coupling fixed point m. d = 1, characterized by the 
exponents in (f2). 

In a recent series of papers (Fogedby ct al.(1995), Fogedby (1998)) we have 
presented a novel approach to the strong coupling features of the noisy Burgers 
equation in one dimension based on a nonpcrturbativc approach in the asymp- 
totic weak noise limit. Since the singular character of the weak noise limit is an 
essential feature of our approach and is already apparent in the linear case for 
A = it is instructive first to consider this case. 



2 The Edwards — Wilkinson Equation 

In the linear case for A = the Burgers equation (1) takes the form of the 
Edwards- Wilkinson (EW) equation (for the slope field) (Edwards and Wilkinson 
(1982)) 

dujdt = vV'^u Vr? , (13) 

i.e., a linear diffusion equation driven by conserved noise. 

Equation (13) is readily analyzed. The time-dependent probability distribu- 
tion for the wavenumber modes Uk = J dxu{x)exp{—ikx), u% = U-k, is given 
by 



P{{uk},t\{ul})^'[[exp 



u 1 \uk — u^le 



(14) 



AL 1 - e-2''fc'* 

where is the initial value for t = 0. For t ^ oo P{{uk},t) approaches the 
stationary distribution (11), using {^/L)Y^^ \uk\'^ = J dxu{x)'^ . 

In a similar way we obtain for the slope correlations in wfc-space, u{kuj) = 
J dxdt exp {iut — ikx) , the Lorentzian line shape characteristic of a conserved 
diffusive mode, 

{u{kw)u{-k - a;)) = ^2 + (^fc2)2 ' (15) 

implying the static correlations {u{x)u{x')) — {A/2v)5{x — x') in conformity 
with (11). In Fig. 3 we have depicted the correlation function. 

Comparing (15) with the scaling form in (7) we also infer the scaling ex- 
ponents ^ = 1/2 and z = 2, characteristic of diffusion and defining the EW 
universality class. Also, noting that the diffusive term in (13) can be derived 
from a free energy F = {1/2) j dxu^, it follows that the EW equation de- 
scribes the fluctuations in an equilibrium system with "temperature" A/2v, 
i.e., Pst = exp [-{2v/A)F]. 

We note already here that the noise strength A plays a special role. Whereas 
A enters linearly in the correlation function (15), the limit A ^ Q appears as 
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Fig. 3. We depict tlic slope correlation function for tlie diffusive mode in the EW case. 
The Lorentzian is centered about a; = and has the hydrodynamical linewidth i>k^ 
vanishing in the long wave length limit. 



an "essential singularity" in the distribution functions P{{uk},t) and -PytX'^) in 
(11) and (14). The special role of the noise as the relevant small parameter is 
more physically recognized by considering the time-dependent correlations 

{u{kt)u{-kt')) = 

[{uu),- {A/2iy)]exp[~{t + t')i'k^] + {A/2u)cxp[-\t-t'\iyk^] , (16) 

(• • ■)i denotes an initial value average. The basic time scale is set by T{k) = l/vk"^ 
which diverges for fc — > 0, characteristic of a hydrodynamical mode. However, 
there is another "time scale" set by the noise, namely a characteristic crossover 
time 

r(Zi)~r(fc)lni . (17) 

For t,t' <C T(fc) the correlations are nonstationary and depend on the initial 
correlation (•••),, whereas for long times t,t' ^ T{k) the correlations enter a 
stationary, time reversal invariant regime and depends on |t — The crossover 
time t(Z\) defines the onset of the stationary regime. For t, t' » r(A;), t(Z\) noise- 
induced fluctuations built up and the system becomes stationary; for Z\ — > 0, 
t{A) oo, and the system never leaves the transient regime, i.e., there is no 
stationary regime. 

Whereas the singular nature of the weak noise limit is easily understood in 
the context of the linear EW equation; it is in fact precisely equivalent to the low 
temperature limit of the Boltzmann factor exp (— F/T), T = A/2i', the situation 
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is more subtle in the Burgers equation where the nonlinear growth term leads 
to a general nonequilibrium system. 

In order to understand the role of the nonlinear term we first consider the 
well-understood properties of the noiseless or deterministic Burgers equation. 



3 The Noiseless Burgers Equation 

The noiseless Burgers equation has the form 

du/dt = vSl'^u + AuVu , (18) 

and describes the transient decay of an interface subject to the damping term 

i^V^7i in combination with the mode coupling or nonlinear growth term \vSlu. 
Correspondingly, for the height field we obtain the deterministic KPZ equation 

f =-V2/^+^(V/.)^ • (19) 

Interestingly, the noiseless Burgers equation (18) can be solved analytically by 
means of the nonlinear Cole-Hopf transformation (Cole (1951) and Hopf (1950)) 



exp 



2v } 



dxu 



(20) 



which maps (18) onto the linear diffusion equation 

— = vy^w , (21) 

and an initial value analysis can be carried out, in particular in the inviscid limit 
— > 0. The basic transient mode structure consists of i) solitons or shock waves, 
connected by ii) ramp modes, and iii) superposed diffusive modes. 



3.1 Linear Modes for A = 

In the linear case for A = the noiseless Burgers equation reduces to the ordinary 
diffusion equation 

— = vW^u , (22) 

which supports decaying diffusive modes 

u{xt) oc exp [—iiJlt + ikx] , (23) 
with gapless quadratic dispersion 

Lol = -iuk^ . (24) 
In Fig. 4 we have depicted the diffusive spectrum 
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Fig. 4. We show the gapless dispersion law for the diffusive modes for the EW case. 
3.2 Soliton and Ramp Modes for A 7^ 

In the presence of the nonhnear mode couphng term the noiseless Burgers equa- 
tion, violating parity symmetry, supports a localized "right hand" soliton mode 
of the sine-Gordon form 

Uk{x) = u^taiihks{x — xq) , (25) 
with center position Xq and amplitude m+. The inverse width 

K = (26) 

depends on the amplitude; in the inviscid limit ^ 0, fc^ ^ 00, and the soliton 
becomes a sharp shock wave. Owing to the Galilean invariance in (8-9) we also 
easily obtain a boosted soliton moving with velocity v. Denoting the boundary 
values at a; = ±00 by u±, we infer the single soliton condition 

u++U- = ~ , (27) 

expressing the soliton velocity in terms of the boundary values u±. The height 
profile corresponding to a soliton has the form 

h{xt) = ^ In cosh [ks{x — xq)] . (28) 
A 

The soliton mode is a stable "dissipative structure" corresponding to determinis- 
tic input at the boundaries (nonvanishing currents) and dissipation at the soliton 
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Fig. 5. We depict a moving "right hand" soliton for the ' noiseless Burgers equation 
together with the associated downward cusp in the height profile. 



center (hot spot). In Fig. 5 we have depicted a moving soliton and the associated 
height profile. 

Furthermore, in the inviscid limit — > (18) supports a negative-slope ramp 
solution 

u{xt) , (29) 

corresponding to a convex parabolic height profile 

Kxt) fx • (30) 
In Fig. 6 we have shown the ramp solution and the associated height profile. 



3.3 Superposed Linear Modes 

The character of the diffusive modes in the presence of a soliton is examined by 
means of a linear stability analysis of Eq. (18). Inserting 

u = Uk + Su , (31) 
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Fig. 6. We show a ramp solution for tlie noiseless Burgers equation together with the 
associated parabolic height profile. 

and solving the associated linear eigenvalue problem we obtain, in addition to a 
localized translation mode, a band of damped diffusive modes 

5u oc exp [—icukt + ikx] , (32) 

with a soliton-induced gap in the dispersion law 

t^fc = -ivik'^ + kl), ks = . (33) 

In Fig. 7 we have shown the diffusive spectrum in the presence of the soliton. 

The soliton modes and associated superposed diffusive modes are thus an 
essential feature of the nonlinear noiseless Burgers equation; this aspect will 
persist in the noisy case. The morphology of the noiseless Burgers equation is 
clear: The equation is basically damped but in the transient regime dissipative 
structures are present consisting of an initial value induced gas of "right hand" 
solitons connected by ramps with superposed linear diffusive modes. With van- 
ishing slope field at the boundaries there is no deterministic current-input and 
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Fig. 7. We depict the diffusive dispersion law in the presence of a soliton for the 
noiseless Burgers equation. The gap is proportional to the soliton amplitude. 



the morphology eventually decays through interaction and coalescence of mov- 
ing solitons. In Fig. 8 we have shown the transient morphology of the noiseless 
Burgers equation and the associated height profile. 




Fig. 8. We show the transient morphology for the noiseless Burgers equation consisting 

of "right hand" solitons connected by ramps for the slope field u and downward cusps 
connected by parabolic segments for the associated height field h. 
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4 The Noisy Burgers Equation 

In Sects. 2 and 3 we discussed the linear EW case in order to illustrate the 
singular character of the weak noise limit and the noiseless Burgers case in order 
to identify the essential nonlinear feature, i.e., the soliton or shock wave, in 
determining the morphology. Here we perform an analysis of the noisy Burgers 
equation (1) including both the nonlinear soliton modes and taking into account 
the singular character of the noise. 

4.1 Phase Space Path Integral Formulation 

The Martin-Siggia-Rosc technique (Martin ct al.(1973)) in a path integral formu- 
lation (Zinn-Justin (1989)) is the appropriate language to analyze the Burgers 
equation and works in the following way: All the information about the statisti- 
cal properties of the Burgers equation is basically contained in the generator or 
characteristic function 



The probability distribution P{{u{x)},t) and correlation functions are easily 
extracted from Z, e.g.. 



Here (• • •) denotes an average withe respect to the noise ij. In order to imple- 
ment the Burgers equation (1) which provides the statistical link between the 
output slope field and the input noise driving the field, we insert the identity 
JY\^^duJS{du/dt — vV'^u — Xu\/u — = 1 in (34)(see note on the Jacobian 
J), expand the delta function, 5{x) oc J dpexp {ipx), and average over the noise 
according to the Gaussian distribution in (2). This procedure yields the path 
integral 




(34) 



{u{xt)u{x't')} = -[5Z 1 5ii{xt)5n{x't')]^=o . 



(35) 




xt 



(36) 



where the action S has the form 




(37) 



with Hamiltonian 



W = - 1 (Vpf + p{uV'^u + XuVu) 



(38) 
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In the linear case for A = we expect the Hamiltonian to reduce to a harmonic 
form. This is accomphshed by performing a shift operation of the noise variable 
p. Setting p = v{iu — ip) we obtain the canonical phase space form 



Z{fj.) oc I ]^rfud(^exp 

•' xt 











i J dxdtjj^u 




exp 




exp 



with bulk action 



and canonical Hamiltonian 



/ 



dxdt 



n 



(39) 



(40) 



(41) 



The shift transformation also generates a surface term Sb which contributes to 
the action for a particular slope configuration, depending on the choice of bound- 
ary condition. It will not play any role in the present discussion; see (Fogedby 
(1998)). 

The Hamiltonian H is composed of a relaxational or irreversible harmonic 
component characterizing the EW case, and a nonlinear reversible mode coupling 
component associated with the drive A in the Burgers-KPZ case. 

The path integral (39) has the usual Feynman form in phase space with u 
and (fi as canonically conjugate variables, u is the canonical "momentum" and 
the noise variable (p the canonical "coordinate" . We note, however, that H does 
not break up in a "kinetic energy" part and a "potential energy" part as is 
characteristic of an ordinary dynamical problem; the coordinate and momen- 
tum variables are mixed in the nonlinear interaction term. Furthermore, the 
Hamiltonian is complex characteristic of the Master equation for a relaxational 
problem. 

Nevertheless, the; Feynman form is very suggestive and allows us to draw 
both on "classical mechanics" in analyzing the path integral and on the associ- 
ated "quantum field theory," yielding the path integral. We note that the noise 
strength A/u plays the role of an effective Planck constant defining a "corre- 
spondence principle" associated with the weak noise limit A ^ 0. 

There are essentially two ways to approach the path integral (39): i) renor- 
malizcid perturbation theory in the nonlinear term A and ii) a nonperturbative 
saddle point calculation for small A. 

Perturbation theory in A yields the same expansion that is obtained by di- 
rectly expanding the noisy Burgers equation (1) and contracting the noise term 
by term according to (3). This expansion which is logarithmically divergent in 
the infrared limit in d = 2 and algebraically divergent below d = 2 requires 
a regularization scheme below d = 2 yielding renormalization group equations 
for I', A, and A. While A is invariant under scaling (yielding the scaling law 
in (10)), the efl[ective expansion coefiicient is X'^A/iy^ yielding the exponents 
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(C, z) = (1/2, 3/2) in d = 1. Note that the expansion is effectively in A and thus 
does not retrieve singular terms in A. 

As discussed earlier the noise A is the small parameter in the problem, at 
least in = 1 where we have carried out the analysis. The singular terms in A 
are precisely obtained by performing a saddle point analysis of the path inte- 
gral, corresponding to the WKB or "quasi-classical" limit of the Feynman path 
integral. 



4.2 Field Equation for A —s^ 

In the weak noise limit A ^ the path integral formulation allows for "a prin- 
ciple of least action" in that the leading saddle point or stationary contributions 
are obtained from the variational principle 6S = with respect to variations 6u 
and d(fi. We thus obtain the Hamiltonian equations of motion 

du 6H ^ dip SH 

where H = J dxH and inserting (41) the field equations 

= -iuV'^if + XuVu (43) 
= +iuV'^u + XuVip . (44) 



du 

dt 
dip 

'dt 



The coupled deterministic field equations (43-44) are an essential result of our 
analysis. In the weak noise limit Zi — > they effectively replace the noisy Burgers 
equation (1). 

Note that we here have obtained a definite level of simplification in handling 

the statistical problem. The stochastic character of the Langevin equation (1) 
has been replaced by the path integral, yielding the additional variable ip charac- 
terizing the noise, and in the weak noise limit the coupled field equations (43-44). 
The alternative formulation in terms of the deterministic Fokker-Planck equa- 
tion is technically more difficult since the Fokker-Planck equation is a functional- 
differential equation, in fact "the Schrodinger equation" for the path integral. In 
Fig. 9 we have depicted the paths in u(p phase space. 



4.3 Linear Diffusive Modes 

It is instructive first to discuss the linear EW case. For A = the field equations 
(43-44) reduce to the linear coupled pair 

— =-zi/VV (45) 
^ = +iyV^u , (46) 
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u', (p' 




u", (p" 



u(p phase space 



Fig. 9. We depict the paths in utp phase space connecting the initial values u' ,Lp' to 
the final values u" , ip" . The fully drawn extremal path corresponds to the stationary or 
saddle point solution for Zi — > 0; the dashed paths characterize the fluctuations about 
the stationary solution. 

admitting the solution 

u{xi) ~ Uq exp [—iuu^t + ikx] + Uq exp [itolt — ikx] (47) 

with gapless dispersion given by (24). We note the doubling of solutions, i.e., 
the equations imply both a damped and a growing solution, unlike the noiseless 
Burgers equation discussed in Sect. 3.1. This is a feature of the noisy case: In 
the stationary regime we attain time reversal invariance and both solutions are 
required in order to describe the stationary correlations. 

For A = the path integral is Gaussian and we leave it as an exercise to derive 
the correlation function (15); it is also not difficult to derive the distribution (14). 

4.4 Nonlinear Soliton Modes 

Leaving aside the question whether the field equations (43-44) in th(^ nonlinear 
case for A 7^ are exactly integrable or admit a Cole-Hopf type transformation, 
we find, like in the case of the noiseless Burgers equation, permanent profile 
soliton solutions, see e.g. (Fogedby (1981)). In the static case we have 

Uo{x) = ±u+t&nhks{x — xq) (48) 

of the same form as (25). Note, however, again a doubling of the solutions like 
in the linear case. In the stationary regime the noise excites both "right hand" 
and "left hand" solitons. Both modes are required in order to correctly describe 
the stationary growth morphology. In Fig. 10 we have depicted the two soliton 
solutions and the associated height field. 



Fig. 10. We depict the "right hand" and "left hand" solitons for the noisy Burgers 
equation and the associated upward and downward cusps in the associated height 
field. 

4.5 Superposed Linear Modes 

Like in the noiseless case we can perform a linear stability analysis of the field 
equations by inserting u = uq + Su and (p = (fo + 5(p. The resulting linear 
eigenvalue problem is exactly soluble, see e.g. (Fogedby et al. (1985)). In addition 
to a translation mode associated with a displacement of the soliton position xq, 
we obtain a band of modes 

Su ~ u'^ exp [—iiVkt + ikx] + u~ exp [iujkt — ikx\ (49) 

with a diffusive dispersion law with a gap given by Eq. (33). In Fig. 11 we have 
depicted a "right hand" soliton with a superposed diffusive mode. 

4.6 A Growing Interface as a Dilute Soliton Gas 

The saddle point modes, i.e., the solitons and diffusive modes, excited by the 
noise immediately allow a physical description of a growing interface. First we 
observe that the saddle point approximation actually implies a dilute gas of 
solitons connected by constant slope segments and satisfying the soliton condi- 
tion (27); in the limit of well-separated solitons these configurations also provide 
saddle point solutions (Zinn-Justin (1989)). A single soliton corresponds to a 
downward or upward cusp in h. A pair of solitons describe a moving step, etc. 
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Fig. 11. We depict a "right hand" soUton (the dashed curve) with a superposed diffu- 
sive mode (the sohd curve); the superposed damped hnear mode is phase-shifted and 
exhibits a gap in the spectrum. 





Fig. 12. We show the stationary growth morphology for the noisy Burgers equation, 
consisting of "right hand" and "left hand" solitons connected by horizontal segments 
for the slope field u and downward and upward cusps connected by constant slope 
segments for the associated height field h. 
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In Fig. 12 we have depicted a general growth morphology in the slope field u 
and the associated height field h. Superposed on the soliton configuration is a 
gas of diffusive modes (not indicated on the figure). We note that the stationary 
growth morphology is quite different firom the transient morphology shown in 
Fig. 8. 



4.7 Dynamics of Solitons 

The Hamiltonian structure of the path integral allows for "a principle of least 
action" and we can associate dynamic and kinetic attributes to the soliton and 
diffusive modes. The energy is inferred from (41), 



E = j dx 



(50) 



The momentum, i.e., the generator of translation, has a form inferred from the 
Poisson bracket {u{x) , (p{x')} = S{x — x') for the canonically conjugate variable 
u and (fi, 



P 



j dxuV(p . (51) 



For a single soliton we have in terms of the boundary values 



E= ±i^[ul-ui] (52) 



P= ±t-[ul-ul] , (53) 



where ± indicates the "right hand" and "left hand" solitons, respectively. 

In particular, for a two-soliton configuration representing a growing step with 
vanishing slope at the boundaries we obtain 

E- 14 (54, 

P = -4:i-^sgn{v) . (55) 
A 

We note the nonlinear velocity dependence characteristic of a soliton mode, see 
e.g. (Fogedby (1981)). Eliminating the velocity we infer the soliton dispersion 
law 

E = iX^\Pf/\ (56) 



and we recover the dynamic exponent z = 3/2 characterizing the gapless soliton 
dispersion. In Fig. 13 we have depicted the soliton dispersion law. 
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33/2 



Fig. 13. We show the gapless sohton dispersion law characterized by the "fractional" 

exponent 3/2. 



4.8 Stationary Distribution 

As alluded to in Sect. 1 it is known that the stationary distribution in the 
Burgers-KPZ case is Gaussian and given by (11) (Huse, Henley, and Fisher 
(1985)). Here wc attempt to make this result, which follows easily from the 
Fokker-Planck equation, plausible within the Martin-Siggia-Rose formalism. 
Prom (36) it follows that the stationary distribution is given by 



JJ dpdu exp 



-S 



(57) 



i{x, 0) and u" = u(x, T) are the initial and final values of the slope 



where u' 

field and S is integrated from t = to t = T. In the weak noise limit A 
variation of S yields saddle point equations equivalent to (43-44), 







— = uV^u + XuVu + iV^p 
at 

^ = -uS/^p + XuVp . 

at 



(58) 
(59) 



Using that the energy vanishes in the stationary state attained in the limit 
T — > oo, we obtain from (37) at the saddle point 



^st(^) ^ 

1 — >-CX) 



\1 



axdtp— 
at 



(60) 
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where p and u are solutions of the field equations (58-59) . In order to demonstrate 
that —ip 2vu for T ^ oo implying the stationary distribution (11), Pg^ oc 
exp [—{v/A) J dxu"^], we define the deviation Au according to —ip = 2i'{u + Au) 
and find to linear order 

Au = i'W'^Au . (61) 

Owing to the Galilean invariancc of the operator d/dt — XuV wc can choose 
an instantaneous frame with vanishing u and (61) implies a decaying solution 
Au oc exp (— z/fc^f). This is not a rigorous proof; basically we assume that the 
trajectories on the E = energy surface in U(p phase space at long times are 
attracted to the sub manifold defined by —ip = 2vu. 

5 Scaling and Universality Classes 

The present path integral approach also permits a simple discussion of the scaling 
properties of the noisy Burgers equation. Focussing on the slope correlation 
function (39) yields 

{u{xt)u{x't')) (X j Y[dud<^exp[i{iy/A)S]u{xt)u{x't') , (62) 

where we assuming vanishing boundary conditions have set Sb — 0. The direct 
evaluation of the path integral (62) requires the application of methods from 
quantum chaos such as periodic orbit theory, see e.g. (Dashen et al. (1974)) 
and is still in progress. However, by discussing (uu) in terms of the underlying 
"quantum field theory" it is an easy task to extract the scaling properties. 

In the "quantum description" the canonical fields u and are replaced by 
"quantum operators" u and ip and likewise the Hamiltonian and momentum 
in (41) and (51). The correlation function (62) can thus be expressed as the 
time-ordered product (Zinn-Justin (1989)) 

(u(xt)u(OO)) cx {0\Tu{xt)u{00)\0) , (63) 

where ii evolves in time according to the "quantum Hamiltonian" (41) and |0) 
denotes the zero-energy stationary state. Displacing the field from {x, t) to (0, 0), 
using the Hamiltonian and momentum operators, we have 

u{xt) = exp [i{u/A){Px + Ht)]u{00) exp [-i{u/A){Px + Ht)] , (64) 

which inserted in (63) together with a complete set of intermediate quasi-particle 
wavcrmmber states \K), P = {A/v)K, with frequency fi, E = {A/u)f2, yields 
the spectral representation 



9 , „ 



{u{xt)u{m)) oc / dKG{K) exp [-i{Qt - Kx)] . (65) 
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Here G{K) is a form factor and f2, K are the frequencies and wavenumbers of 
the quasi-particles in the theory. 

The scaling Hmit for large x and large t corresponds to the bottom of the 

quasi-particlc spectrum and wc note that only gapless excitations contribute. 
Assuming a general dispersion law with exponent /3, 

n = AK^^ , (66) 

and assuming that the form factor G{K) is regular for small wavenumber, 
G{K) cx const., we obtain, rescahng K 

{u{xt)u{00)} oc x-^ J div:e-i^^^(*/^^)-i^ . (67) 

Comparing the spectral form (67) with the dynamic scaling hypothesis (7) we 
first infer the robust roughness exponent C = 1/2, independent of the quasi- 
particle dispersion law. The dynamic exponent z is given by the exponent /3 
for the quasi-particle dispersion law. In the linear EW case the gapless diffusive 
dispersion law f2 oc yields the dynamic exponent 2; = 2; in the Burgers-KPZ 
case the diffusive modes develop a gap and do not contribute to the scaling; 
however, the noise excites a new nonlinear gapless soliton mode with dispersion 
i7 oc K^/"^, yielding the exponent z = 3/2. We also obtain a heuristic expression 
for the scaling function F3 in (7), 




which has the same form as the probability distribution for Levy flights (Fogedby 
(1994)). In Table 1 we have summarized the exponents and universality classes 
for the EW and Burgers-KPZ cases. 



Table 1. Scaling exponents and universality classes 



Model 


Roughness exp. {() 


Dynamic exp. (z) 


Universality class 


EW equation 


1 
2 


2 


EW 


Burgers equation 


1 
2 


3 
2 


Burgers-KPZ 



6 Summary and conclusion 

We have here advanced a novel approach to the growth morphology and scaling 
behavior of the noisy Burgers equation in one dimension. Using the Martin- 
Siggia-Rose (MSR) technique in a canonical form we have demonstrated that 
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the physics of the strong coupHng fixed point is associated with an essential 
singularity in the noise strength and can be accessed by appropriate theoretical 
soliton techniques. 

The canonical representation of the MSR functional integral in terms of a 
Feynman phase space path integral with a complex Hamiltonian identifies the 
noise strength as the relevant small nonperturbative parameter and allows for a 
principle of least action. In the asymptotic weak noise limit the leading contri- 
butions to the path integral are given by a dilute gas of solitons with superposed 
linear diffusive modes. The canonical variables are the local slope of the interface 
and an associated "conjugate" noise field, characteristic of the MSR formalism. 
In terms of the local slopc^ the soliton and diffusive mode picture provides a many- 
body description of a growing interface governed by the noisy Burgers equation. 
The noise-induced slope fluctuations are here represented by the various paths or 
configurations contributing to the path integral. Moreover, a spectral represen- 
tation of the slope correlations based on the underlying "quantum field theory" 
gives access to the scaling exponents and provide a heuristic expression for the 
scaling function. 

So far the present nonperturbative approach has only been implemented in 
the one dimensional case where the analysis is tractable. However, the singular 
nature of the weak noise limit seems to be a general feature of the onset of the 
stationary regime and might also be important for the Burgers-KPZ equations 
in higher dimension. It remains to be investigated whether the present saddle 
point approach can be generalized to this case. 
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